The electron effective mass in the strongly correlated 2D-uniform electron fluid from 

finite-temperature calculations. 
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The very-low temperature thermal effective mass m* of paramagnetic and ferromagnetic electrons 
in a uniform electron fluid in two dimensions is studied. Analytical and numerical evaluations are 
used to meaningfully define an m* even in the the Hartree-Fock approximation. The Hartree-Fock 
m* decreases linearly with the electron-disk radius r s . Correlation effects lead to strong cancellations 
between exchange and correlation. Thus the effective mass is enhanced with increasing r s for the 
unpolarized fluid, while m* decreases with the r a of the polarized fluid. The effective mass is 
calculated from the coefficient of the quadratic temperature dependence of exchange-correlation 
free energy F xc . This is calculated in a physically transparent manner using a new formula for 
the effective mass. This uses the T — pair-distribution functions of Gori-Giorgi et al., and the 
temperature derivative of a quantum analogue of the potential of mean-force well known in the 
statistical mechanics of classical fluids. The results are compared with recent quantum Monte-Carlo 
simulations at T = 0, as well as with other available experimental and theoretical data for the 
effective mass. 

PACS numbers: PACS Numbers: 05.30.Fk, 71.10,+x, 71.45.Gm 
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I. INTRODUCTION. 

At very low temperatures, essentially at T = 0, the 
electrons in a 2D uniform electron fluid (UEF) occupy 
momentum states up to the Fermi energy Ep. If the elec- 
trons were non-interacting they have only kinetic energy, 
the potential energy being a constant. The energy bands 
of the electrons are of the form k 2 /2m, where the mass 
m of the non-interacting electrons is unity, even in semi- 
conductor nanostructures or metals where the material 
parameters can be subsumed by a suitable re-definition 
of the effective atomic units [if. Such 2D electron lay- 
ers have been experimentally realized in semiconductor 
interfaces, enabling the experimental probing of such sys- 
tems which are of great fundamental and technical im- 
portance. 

If the the number of electrons in a unit area (mea- 
sured in atomic units) is n, the radius of the disk in real 
space per electron, viz. r s , is l/y/(wn). The ratio of the 
Coulomb interaction to the kinetic energy is a measure 
of the 'strength' of the interactions in the system. In the 
limit T — t it is found to be equal to r s . Thus r s plays 
the role of the 'expansion parameter' or coupling con- 
stant in the perturbation theories of the electron fluid. 
Perturbation theories work reliably when the coupling r s 
is less than unity. The range of validity can sometimes 
be increased by the use of a 'screened interactions', at 
least for some properties. 

One of the tenets of Fermi liquid theory, originally 
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developed to treat the behaviour of electrons in metals 
(where usually 2 < r s < 5 ), is that the low-energy ex- 
citations in these interacting systems, known as quasi- 
particles, are very similar to non-interacting electrons. 
However, the mass m is replaced by an effective mass m* 
which allows for the modification of the energy bands 
due to many-body interactions. Other parameters which 
describe the electrons, e.g., the Lande g factor, are also 
modified and these are the 'Fermi-liquid parameters' of 
Landau theory. Unfortunately, standard many-body the- 
ory which uses perturbation methods cannot make a reli- 
able evaluation of m* for systems where r s is greater than 
unity. Physically motivated approximations are needed 
to truncated the perturbation chain, and these invari- 
ably result in the failure to satisfy the sum rules, self- 
consistency conditions etc. If the results were insensitive 
to the various possible choices for screening, vertex cor- 
rections, selection of graphs, etc., then this would not 
matter. Unfortunately, the evaluation of m* , usually 
carried out from a physically motivated approximation 
to the self energy, turns out to be very sensitive to the 
model used. Thus, for example, perturbation theories 
may predict the m* of electrons at r s =5 to be a fraction 
of the bare mass, or several times the bare mass, depend- 
ing on the model used. These methods even fail to predict 
positive-definite pair-distribution functions (PDFs), i.e., 
g(r) for useful values of r s . 

However, significantly more reliable results are avail- 
able from quantum Monte Carlo simulations (QMC), 
especially for PDFs and correlation energies. Unfortu- 
nately, the m* evaluations involve the probing of excited 
states, and only a few QMC calculations are available for 
the effective mass of the 2D UEF. Experimentally too, 
the measurements have been very challenging. Neverthe- 
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less, currently available results, both experimental and 
simulational, show that the to* of the 2D-paramagnetic 
electron liquid is enhanced above unity as r s increases, 
while the to* of the fully spin-polarized system remains 
below unity even as r s is increased. 

An alternative approach to the study of the effective 
mass is to look at quasiparticle excitations in a fluid at a 
temperature T close to T = 0. The thermal excitations 
occur in a strip of energy of width T near Ep. Th ex- 
citations in the system are associated with an increase 
in the Hclmholtz free energy Fq (non-interacting case), 
or F (interacting case), of the system. The quasiparticle 
mass enters directly into the specific heat of the system. 
Hence the ratio of the specific heats of the interacting 
system and the non-interacting system provides a direct 
and unambiguous measure of the effective mass m* of the 
excitations. Under certain conditions, these can be iden- 
tified with the Landau quasiparticles of the Fermi liquid 
in the limit T — > as discussed by Luttinger and others 
Thus a calculation of the interacting free energy F 
as a function of r s at finite-T would provide an estimate 
of the effective mass to* , if the finite-T free-energy of the 
system could be calculated. 

The calculation of the finite-T interacting free energy 
or the self-energy at the Hartree-Fock (HF) level is well 
controlled. However, the self-energy is divergent at the 
Fermi energy and it is not possible to define an effective 
mass. On the other hand, the logarithmic divergencies 
in the HF-free energy can be separated out and we are 
able to present a thermal effective mass at the HF level, 
viz., to*, arising entirely from exchange processes. Here 
we present previously unpublished results for the coeffi- 
cients of the T 2 term of the polarized system. The HF to* 
decreases with r s and this approximation breaks down for 
r s > 3 when m x becomes negative. Thus the inclusion 
of correlation corrections is imperative to obtain a phys- 
ically meaningful result. 

Instead of using diagrammatic methods, in our previ- 
ous work we used a calculation of F xc at finite-T via a 
coupling-constant integration of the PDFs of the interact- 
ing system. The finite-T g{r) needed for the calculations 
were obtained using the classical-map hyper-netted-chain 
(CHNC) method where the 2D electron system at the 
temperature T = is replaced by a classical Coulomb 
fluid at the temperature T q . Given that the interacting 
chemical potential /i becomes negative beyond r s ~ 2, 
(e.g., at r s = 5, [i/Ef ~ — 5), the distribution func- 
tions even at T just slightly above T = are classical 
Boltzmann distributions with little or no occupation at 
ftp. Hence the study of a classical model which correctly 
incorporates quantum features via effective potentials is 
a very reasonable proposition. However, although the 
CHNC can be unambiguously implemented at T = and 
at sufficiently elevated temperatures, we are beset with a 
number of difficulties in dealing with the 'warm-dense' re- 
gion close to T = 0. The two main difficulties are (a) the 
elimination of logarithmic-divergent terms which need to 
exactly cancel with the exchange and correlation contri- 



butions, (b) the finite-T modeling of the bridge function 
of the classical-fluid which controls the cluster diagrams 
beyond the hyper-netted-chain sum of diagrams. In this 
study we present an alternative approach which partially 
circumvents these difficulties, and provides a more trans- 
parent analysis, leading to a new formulation of the ef- 
fective mass calculation. 

The g(r s ,r), i.e., PDFs of the 2D-electron system (at 
T = 0) have been accurately parametrized by Giri-Giorgi 
et al., Q and may be considered known. The F xc at 
T = can be written as a coupling-constant integration 
over g(Xr s ,r) where A is the coupling constant. Further 
more, we consider -log{g(r)} as the potential of mean 
force, viz., /3V m f(r) of the equivalent classical interacting 
Coulomb fluid. Then we determine the finite-T form of 
/3V m f (r) to second order in the temperature and use this 
to directly evaluate the second-order temperature cor- 
rection to the exchange-correlation free energy. In the 
following we show that the method leads to a transpar- 
ent, if approximate, calculation of the effective mass m* 
at arbitrary polarizations, and in good agreement with 
the available results in the field. 



II. THEORY 

The thermal effective mass to* can be expressed as a 
ratio of the heat capacities of the interacting and non- 
interacting systems as: 



[d 2 F{T)/dT 2 ' 
[d 2 F a {T)/dT 2 



(1) 



Here C v is the specific heat at constant volume. The in- 
teracting free energy F is the sum Fq + F x + F c , where 
F x and F c are the exchange and correlation contributions, 
with F xc = F x + F c . Hence the problem of determining 
to* reduces to a calculation of exchange-correlation ef- 
fects at finite temperatures, near T = 0. However, such 
calculations are in many ways even more demanding than 
those at zero temperature, as perturbation methods have 
to now deal with a whole host of new diagrams, their sin- 
gularities and cancellations. These difficulties were first 
addressed in the papers by Luttinger, Ward, and Kohn 
2]. Similarly, QMC methods are also equally difficult, 
especially for T very close to zero. 

The Hartree-Fock self-energy becomes logarithmically 
divergent near kp and it is not possible to define an to* 
via the self-energy at T = 0. In the following we first 
examine the exchange-only, i.e., Hartree-Fock, approxi- 
mation to F at finite temperature, and calculate a regu- 
larized effective mass to* = tyihf = 1 + Ato^ for para- 
magnetic and ferromagnetic 2D electrons, i.e., for spin 
polarizations £ — and 1. It is found that m x , con- 
taining the corrections from F x , is a decreasing linear 
function of r s for both polarizations. 
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A. Non-interacting and Hartree-Fock Helmholtz 
free energies 

The non-interacting free energy Fq = Eq — TS, where 
S is the entropy, can be evaluated and expressed in terms 
of the dilogarithm function. Writing 77 = fx/T, where fi is 
the chemical potential, Fq and Eq for the 2D system are 
given as follows, per unit atomic volume, and in Hartree 
atomic units: 
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Alternative expressions are possible, as given in ref [3, Q ■ 
The low-temperature expansion of F (t) can be given as: 

F (t,C = 0) - £ (0,0)(l- 3.2859i 2 ) C = (4) 
F (t,C=l) - S (0,l)(l-0.82246t 2 ) C = l 



Eo(P,Q/n = 0.5(1 + C)/r 2 s 



(5) 
(6) 



The quadratic coefficient in the expansion of Fo(t, C) 
would be denoted by 0.2(C) where needed. Thus a 2 (C = 
1) = — 0.82246i?o- Note that here we have used the same 
paramagnetic Ep in defining t = T/Ef for all 2D sys- 
tems mentioned in the above equations. 

The first-order (i.e., unscreened) exchange free energy 
F x consists of F£, where i runs over the two spin species. 
At T = these reduce to the exchange energies: 



Ef/n = 



30r 



1/2 



(7) 



Here m = n(l + Q/2, and n 2 = n(l - Q/2. Then the 
exchange energy per particle at T = 0, i.e., the internal 
energy contribution E x /n at T = becomes 



E x /n={Ef+E%)/n 



3/ 2 , 3/ 2 l 



37rr s 



(8) 



Here c\ and C2 are the fractional compositions (1 ±0/2 
of the two spin species. 

We also define the species-dependent reduced chemical 
potentials fi®/T by rji, reduced temperatures t\ = t/(l + 
Q and £2 = t/(l — Q, based on the two Fermi energies 
Epi and Ep2 which are Ep(l ± (). Then we have: 



w * " 16 1 loo {^-uyn 

The I-1/2 is the Fermi integral defined as usual: 



da; a;" 
'0 l + e 31 " 



The r/i are given by 



77, = logCe 1 /'* - 1) 



(9) 



(10) 



(11) 



In the paramagnetic case Eq.[9]reduces to the result given 
by Isihara et al. [(| (see their Eqs. 3.4-3.6; they use a 
slightly different definition of the Fermi integral). 

The total exchange free energy per unit atomic vol- 
ume is F x = YiF? . The accurate numerical evaluation of 
Eq. [3] requires the removal of the square-root singularity 
by adding and subtracting, e.g., I 2 (— \t]\)/(v — H) 1 / 2 for 
the case where rj is negative, and v = u, and so on. 

A real-space formulation of F x = Ff + F'£ using the 
zcroth-order PDFs fits naturally with the approach of 
our study using pair-distribution functions of the electron 
fluid as the main ingredient. Thus 



(12) 



Here h^(r) — g%(r) — 1. In the non-interacting system 
at temperature T, the antiparallel hi 2 , viz., gi 2 (r, T) — 1, 
is zero while 
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/i?i(r) = ^T, klik2 n(ki)n(k 2 )e 



:(ki-k 2 )T = _ 



= ~[f(r)Y 



Here k, r are 2-D vectors and n(k) is the Fermi oc- 
cupation number at the temperature T. At T = 
/(r) = 2Ji(kir)/k r where J\{x) is a Bessel function. As 
a numerical check, we have evaluated the exchange free 
energy by both methods, i.e., via fc-space and r-space cal- 
culations. 

The following small-T expansions are useful for our 
purposes: 

F x (r s ,t,( = Q) = E x (r s ,C = 0)[l + (7T 2 /16)t 2 log(t) 

••] (13) 



-0.56736^ 



This result, i.e., F x for the unpolarized system, has been 
given by Isihara et al., Q and re-confirmed by Mahan et 
al ■ The corresponding internal energy, E x can be ob- 
tained from the relation E = d{/3F}/d(3 where (3 = 1/T. 
We refer to the coefficient of the t 2 log(t) terms as Al, 
and that of the quadratic term as A 2 . These depend 
on the spin polarization £. The results for the fully po- 
larised case do not seem to be previously available in the 
literature. We have obtained the following expansion: 

F x (r s ,t,( = l) = E x (r s X=l)[l + (7T 2 /64)t 2 log(t) 

(14) 



-0.28368t 2 



B. The Hartree-Fock effective mass 

The Hartree-Fock electron fluid is a theoretical con- 
struct which does not exist in nature. However, it is a 
very useful conceptual model. The difficulties in this con- 
ceptual model are due to the singular behaviour of the 
single-particle energy, self-energy etc., close to the Fermi 
energy. Thus taking derivatives near the Fermi energy 
becomes meaningless. However, as these divergences are 
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M m x (?=0) 
B-B m* x (^=l) 
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FIG. 1: The panel (a) shows the modification of the chemi- 
cal potential /i due to interactions, driving it to the negative 
regime for r s >~ 1. Negative fi is typical of classical fluids, 
(b) The Hartree-Fock effective mass, m*, is shown as a func- 
tion of the density parameter r 3 , for the polarizations £ = 
and 1. 



logarithmic, they can be integrated over and the free en- 
ergy and related quantities can be evaluated. They are 
found to contain logarithmic terms which are explicitly 
exposed in Eqs. Q2] and Q31 These logarithmic terms, of 
the form A L (()t 2 log(t) are removed when higher order 
corrections are included in the theory, while the quadratic 
term A2(()t 2 contributes to the heat capacity. Hence it 
is clear that we can define a regularized effective mass 
by dropping the t 2 log(t) terms and calculating an m x 
entirely from the coefficient A 2 . Thus 



m x (0 = l + A 2 (0/a 2 (0 



(15) 



Here a 2 is the coefficient defined in Eq.0]in the expansion 
of Fq. Results from such calculations are shown in Fig.Q] 
These results show how the exchange contributions widen 
the Hartree bandwidth with the decrease of the effective 
mass. For r s > ~ 3, the value of m* becomes negative. 
Clearly, the Hartree-Fock model is misleading except at 
smaller r„. 



III. CORRELATION CORRECTIONS TO THE 
FREE ENERGY AND THE EFFECTIVE MASS. 

Microscopic theories of Landau Fermi liquids calculate 
m* from the solutions of the Dyson equation for the one- 
particle interacting Green's function of the system. If the 
real part of the retarded self-energy is £i(fc, w), the Lan- 
dau quasi-particle excitation energy Eqp(Ic), measured 



with respect to the chemical potential is used in calcu- 
lating the effective mass to*. 



E QP (k) 
1 

to* 



e k + Ei(fc,w)| tJ= _E QP 
k 2 /2 - E F 
dE QP (k) 



kp dk 



\k=k f 



(16) 
(17) 

(18) 



The quasiparticle energy is the eigenvalue of the Dyson 
equation and contains the non-interacting energy e{k) 
plus the real part of the the selfenergy E(fc,w) evalu- 
ated self-consistently at the quasiparticle energy itself. 
We may rewrite the selfenergy as an exchange part, and 
a correlation contribution: 



E = E i; (fc) + E c (A;,£ 



(19) 



This should be evaluated self-consistently to satisfy sum 
rules, Ward identities etc., to give a conserving approx- 
imation. In practice, these self-consistency conditions 
have to be stringently satisfied if the logarithmic singular- 
ities in Yi x are to be completely cancelled by correspond- 
ing singular terms contained in E c . Thus any approxi- 
mate theory must necessarily be on guard against spuri- 
ous contamination from inadequate cancellations. These 
could greatly enhance the calculated to*, as is indeed 
observed in some RPA-type calculations. Further more, 
the interacting chemical potential fJ,/Ep is very different 
from unity, and strongly negative. These make the at- 
tempt to impose self-consistency a very daunting task. In 
the temperature theory, say Eq. [TS] of the previous sub- 
section, these would manifest as spurious contributions 
to m* proportional to l/t as t — > 0. 

The exchange and correlation free energy F xc (r s ,t) 
at any temperature and spin-polarization can be eval- 
uated from pair-distributions functions rather than from 
self-energies or Green's functions. The exchange free 
energy F x is just the contribution to a coupling con- 
stant integration (see below) over the PDFs at zero cou- 
pling (see Eq. IT2"j) . The F x , evaluated from the non- 
interaction PDF, g°(r) via Eq. [T^] contains singular loga- 
rithmic terms. These singular terms at zero coupling are 
offset by the contributions from the rest of the coupling- 
constant integration reaching out to full coupling. 



F X c(rs,t)/n 



dX- 



2nrdr 



E 



CiCj(gij(r, A) - 1) 



(20) 

Here A is the coupling constant. The r s , and t = T/Ep 
dependencies in the (A, r) are not displaced for brevity. 
In our previous work (e.g., ref. [8| and references there-in) 
the needed gtj (r) are calculated at any given temperature 
T (including T = 0) using the CHNC method for 2D sys- 
tems [4, 5;]. In CHNC, the quantum fluid at T = is re- 
placed by a classical fluid at T q = 1//3. At finite temper- 
atures T > 0, f3 = l/^jT 2 + T 2 . Although this method 
gave reasonable results for m* via Eq. Q] the question 
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of the extension of model bridge functions to finite-T, 
as well as the accuracy of the cancellation of logarithmic 
terms which should cancel accurately remained troubling 
issues. Numerical calculations very close to T = are 
also very susceptible to difficulties due to the sharpness 
of the Fermi functions at very low T. 

Accurate pair-distribution functions of the 2D-UEF at 
T = are now available from QMC as well as from ana- 
lytical representations developed by Gori-Giorgi et al 
We use these functions instead of those generated from 
the CHNC, as this is equivalent to replacing the hard- 
disk bridge functions of the CHNC description with the 
Coulomb bridge function [9|. If the 2D-PDF, i.e., g(r), 
obtained by these methods were that of a classical fluid, 
then it would be of the form: 



j(r) H> exp [—/3V m f(r)] 



(21) 



V m f(r) is known as the potential of mean force, and 
is simply the Kohn-Sham potential at r in a classical 
fluid where one particle is already at the origin. The 
CHNC attempts to construct f3V m f(r) directly from the 
diffraction-corrected Coulomb potential Vd c (r) and the 
Pauli-exclusion potential P(r)S ss i via the modified HNC 
equation. These potentials are discussed in greater detail 
in refs. 0, @. 

g(r) = exp [-fiV dc (r) - (3P{r) + N(r) + B(r)} (22) 

Here N(r) is known as the nodal function, while B(r) 
is the bridge functions [10| which bring in multi-particle 
clustering effects which are outside the scope of the 
hyper-netted-chain diagrams. The potentials Vd c (r), 
and P{r) are long-ranged, while N(r),B(r) contain 
many-body effects, and screening effects which damp 
Vdc(r), P(r), so that g(r) —> 1 for large r. 

In the present paper we follow a different strategy. The 
potential of mean force, j3V m f(r,T) at finite-T is devel- 
oped as a Taylor expansion around T = 0, i.e., around 
(3 = 1/T q and F xc is calculated using the T — PDF. 



P 

PV m} {r,T) 

F xc {T)/n 
AF~ XC 



IKT^ + T) 1 / 2 
/3V m/ (r,0) + (T72) 



(23) 



2^d 2 pV mf (r,T). 



dT 2 



T=0 



(24) 



F xc {T = 0)/n+(T 2 /2)AF xc 



(TV2)/dA 2/ 

d 2 (3V mf (r,T, A) 



r(r,T,A) (25) 



T(r,TX) = ^2ciCjgij(r,X) 



dT 2 



\T=0 



In calculating the second-order temperature derivative 
of the potential of mean force, we use the CHNC form, 
Eq. I22[ and construct an approximate simplified from as 
a screened pair potential. 



The effect of the nodal term N(r) and the bridge term 
of Eq. [55] is simplified and subsumed in the screening 
wavevector k sc . At high densities, i.e., r s < 1, elec- 
trostatic potentials screen at the Thomas-Fermi value 
2irn/Ep = 2, and drop to about half this value at low 
densities. This is further discussed in the section on re- 
sults, given below. Also, the T-dependence or polariza- 
tion dependence in the screening parameter is neglected 
in the present study. 

The diffraction corrected Coulomb potential and the 
Pauli exclusion potential are of the form 



V dc = (l/r){l - e- k *» r } 



(27) 



The de Broglie thermal wave vector kth = { 7r T;} 1 ^ 2 at 
T = appears in Vd c (r). Its temperature dependence 
is neglected in our context as it only modifies the very 
small-r region, where g(r) is negligible. Hence it plays no 
significant role in Eq. 1251 The Pauli exclusion potential 
j3P{r) is given by 



pP(r)S s 



-lo 9 {g°(r)} + (g°(r) 



1 



c °(r) (28) 



pV mf (r,T) ~ {l3V dc (r) + /3P(r,T)}e k ^ 



(26) 



Here c°(r) denotes the 'direct-correlation function' of 
Ornstein-Zernike theory. It can be calculated from g°(r), 
i.e., the non-interacting pair-distribution function for 
parallel spins s — s'). An average potential may also 
be defined for use with the paramagnetic system treated 
as a single species, using g® v = 0.5(^1 + .912) m Eq. [25] 
Hence P(r) is known at any temperature and r s . It is a 
universal function of r/r s , and known only as the product 
of j3 and P(r). However, we may use the CHNC values 
of T q to define a (3 at any r s and determine P{r) by 
itself. We have numerically calculated its second-order 
temperature derivative (Fig. [2]) and denote it as f3P( 2 ) 
in the following. These quantities are now assembled to- 
gether to calculate the second-order finite-T correction 
to F XC (Q. We have: 

T(r,TA) = T(mf) + T(P) (29) 
T(mf) = ^cc,^(r,A) ^ r ' 0) ( 30 ) 

ij J 9 

r(P) = -J2c l c m (r,X)l3P {2 Hr,t)S l] \ t=Q (31) 

ij 

m*(C) = l + AF xc (0/a 2 (0 (32) 

The 02(C) coefficient in the last equation is the coefficient 
of the T 2 term in the interacting free energy, as in Eq. 0] 
The first term in the r.h.s. of Eq. [29] arises from the 
temperature derivative of f3 in /3V m f (r, t) , evaluated at 
T = 0, while the second term is from the second temper- 
ature derivative of the Pauli exclusion potential. In eval- 
uating Eq. [22 we use the pair-distribution functions of 
Giri-Giorgi et al. [3], as the approximate forms derived 
from CHNC using a hard-disk bridge function are less 
accurate. The 2D quantum temperature T q and other 
CHNC procedures used here are those given by Perrot 
and Dharma-wardana in ref . Q . 
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FIG. 2: The Pauli exclusion potential f3P[r) for 2D-electrons, 
and its second temperature derivative (t = T/Ef), derived 
from <7ii(r). The bottom panel shows the results for the 
paramagnetic, i.e., averaged g°(r). These potentials exactly 
reproduce the Fermi hole in the pair-distribution functions. 



Eq. [21] provides a new model for the calculation of 
the effective mass to* as a sum of contributions from 
the temperature dependence of the potential of mean 
force, m xc (mf), and the temperature dependence of the 
Pauli-exclusion potential, m xc (P). These contributions 
are shown in Pigs. [3] and H] and labeled according to the 
following equation. 



rn* = 1 + m xc = 1 + m xc (mf) + m xc (P) (33) 
m* = m{mj,P); m* (mf) = 1 + m xc {mf) 



The contribution m xc (mf) may be thought of as re- 
sulting from temperature dependent modifications in 
exchange-correlation contributions as a function of r s 
due to screening, while m(P) are contributions from the 
screened exchange interactions. Only the second-order 
temperature derivatives have been used as the first-order 
derivatives (which contain logarithmic terms of the form 
t 2 log(t) are considered to cancel themselves out in an 
accurate theory. 
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FIG. 3: The effective mass m* as a function of the density 
parameter r s , for unpolarized (£ = 0) 2D electrons. See text 
for more details. 



IV. RESULTS AND DISCUSSION 

The values of to* calculated using different diagram- 
matic perturbation expansions (llL fl2j ] differ significantly 
from each other, and from QMC results [F3T - fl6j ]. How- 
ever, the general trend seems to be that to* increases 
with r s for unpolarized (( = 0) electrons, while the m* 
for C — 1 decreases with r s . We have summarized some 
representative results as well as our results in Fig. [3] for 
C = 0, and Fig. [4] for fully polarized electrons. Asgari 
et al. have given calculations where the Dyson equation 
is solved "self-consistently" (GW-SC), as well as where 
an "on-shell" approximation (OSA) has been used 
They have included approximate vertex corrections us- 
ing the Kukkonenen-Overhauser (KO) approach, but it 
is not clear if the fully interacting chemical potentials (see 
Fig.[T]) have been used. We show the results for to* (to/), 
as well as for the full to* which includes the contribution 
from the screened Pauli exclusion term. In the case of 
( = we have used the temperature derivative of the 
averaged Pauli potentials (lower panel, Fig. [3J where as 
the derivative of the full Pauli potential was used with 
polarized electrons. Te temperature derivates were cal- 
culated from numerical data in the range T/Ep = 0.05 
to 0.15. 

The essential weak point of this calculation is in the 
specifying of the screening factor k sc in Eq. [2H1 A one- 
parameter screening approach is indeed very approxi- 
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FIG. 4: The effective mass m* as a function of the density 
parameter r a , for polarized (£ = 1) 2D electrons. See text for 
more details. 



mate. However, this is guided by our calculation of 
m*(mf) shown in Fig.s[3]and|U While the high-density 



Thomas-Fermi k sc is 2 for electrostatic interactions, we 
have used 2.5 for r s below unity for the screening of the 
exchange interactions brought in by the Pauli potential. 
The screening parameter k sc drops to 2 at r s — 4, then to 
1.54 at r s — 10, and unity near r s = 30. By r s — 30, m* 
becomes 1.34 for C— 0, and takes the value of 0.817 for 
£ = 1. Drummond and Needs have recently reported a 
QMC value of m* = 1.34 at r s = 10 for C = 0. Our model 
does not exclude such higher values of m* , but this would 
require a significantly weaker screening of the potentials, 
which accomodates less comfortably with the behaviour 
of the potential of mean force (eq. [5T] as a function of r s . 



V. CONCLUSION 

We have presented results for the effective mass m* as a 
function of the density and polarization of a 2-D electron 
fluid, using an entirely new approach based on the known 
pair-distribution functions of the 2D fluid. Good agree- 
ment with the trends and magnitudes observed in recent 
experimental and theoretical approaches to the problem 
are recovered. The essential inputs to the present calcu- 
lation are: (i) the pair distribution functions of the 2D- 
electron fluid as parametrized by Gori-Giorgi et al [3[, 
(ii) the Pauli exclusion potential extracted from the non- 
interacting PDF, (iii) elementary thermal physics, and 
(iv) some classical-map HNC concepts from the theory 
of classical statistical mechanics. The method provides a 
physical understanding of the processes contributing to 
the enhancement or diminution of the effective mass. 
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